Enhancing robustness and immunization in geographical networks 
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We find that different geographical structures of networks lead to varied percolation thresholds, 
although these networks may have similar abstract topological structures. Thus, the strategies for 
enhancing robustness and immunization of a geographical network are proposed. Using the gener- 
ating function formalism, we obtain the explicit form of the percolation threshold q c for networks 
containing arbitrary order cycles. For 3-cycles, the dependence of q c on the clustering coefficients is 
ascertained. The analysis substantiates the validity of the strategies with an analytical evidence. 
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Complex networks (see reviews [1]) provide powerful 
tools to investigate complex systems in nature and soci- 
ety. The properties of complex systems are affected by 
the geographical distribution of the components. For ex- 
ample, routers of the Internet Q and transport networks 
3 lay on the two-dimensional surface of the globe; world- 
wide airport network is confined by the geography 
neuronal networks in brains [5( occupy three-dimensional 
space. Thus it is helpful to study the geographical com- 
plex networks [1, 0, S, i, [M El • 

From the abstract geometrical point of view, an ab- 
stract set can describe a general system. When it is 
equipped with some geometrical structures, the set can 
further describe a specified system. So, an abstract topo- 
logical network is an abstract set that consists of nodes 
and links. A metric can be added to the abstract topo- 
logical network, the metric could be arbitrary, not neces- 
sarily the Euclidean metric [6]. Embedding a particular 
abstract topological network, into a suitable metric space 
provides a method to add a metric. 

In this paper, the problem of percolation thresh- 
olds in geographical networks is studied. Three types 
of geographical networks are investigated: the normal 
model 0, @], the hollow model and the concentrated 
model. By extensive numerical simulations, we found 
that the percolation threshold q c (the point that a span- 
ning connected cluster emerges) for these models sat- 
isfies: q c (concentrated) >q c (normal) >q c (hollow). Based 
on these results, we suggest a strategy (hollowing) for 
enhancing robustness and a strategy (concentrating) for 
immunization of geographical networks. The different ge- 
ographical networks have different distribution of cycles. 
The geographical dependence of the percolation thresh- 
old is investigated by the generating function process in 
abstract networks containing cycles. 

Based on the lattice embedded model Q , the networks 
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are generated as follows: each node in an L x L lattice 
with periodic boundary conditions is assigned a degree 
quota k, drawn from the prescribed degree distribution. 
Here, scale-free (P(k) ~ fc~ A ,m ^ k) and exponential 
(P(fc) - e- fc / fc °,m < k) de gree distributions are consid- 
ered, where m is the minimum degree a node can have, 
(a) The normal case [lattice embedded scale-free (LESF) 
or lattice embedded exponential (LEE)]: A node i con- 
nects to its closest neighbors until its degree ki is realized, 
or up to a cutoff distance A^/kl , where A is large enough 
to ensure that almost all the degree quotas can be ful- 
filled, (b) The hollow case [hollow LESF (HLESF) or 
hollow LEE (HLEE)]: similar to the normal case, except 
that a node i has probability p to be forbidden to connect 
its first n nearest neighbors, (c) The concentrated case: 
A is set smaller than that in the normal case. The pro- 
cess is repeated throughout all the nodes in the lattice. 
For the hollow case, when the network degenerates to a 
lattice and n is small, it becomes similar to the tunneling 
effect on Euclidean lattices [l2[ • In the following simula- 
tions, we choose A = 7, n — 8, and p = 1; network size 
N = 10 6 , minimum degree m is 4 for scale-free networks 
and 1 for exponential networks, and all the data are av- 
eraged over 1000 ensembles, unless otherwise specified. 

The algorithms of Newman and Ziff [l3T | is performed 
to calculate the threshold q c , which is defined as the point 
where the differential of the size of the largest cluster as 
a function of occupying probability q maximizes. The 
definition is equivalent to the usual definition by the 
emergence of a spanning cluster for large network size 
limit [14]. Figure [1] shows a clear drop of the perco- 
lation threshold q c in the hollow networks than in the 
normal networks. So the robustness can be enhanced 
by the model. The size effect (scale-free degree distribu- 
tion) is demonstrated in Fig. [5] The LESF and even the 
HLESF networks also have non-zero percolation thresh- 
olds for A S (2,3), the results are consisitent with Ref. 
[l4j. Again, the drop in q c for hollow networks is appar- 
ent. 

The amendment in hollow networks is small in the 
physical space. As Fig. [3] shows, the average spatial 
length (I) of the edges for the hollow networks dose not 
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FIG. 1: (Color online) Percolation thresholds of the networks 
for both random failures (left panels) and intentional attacks 
(right panels). In each subgraph squares represent for nor- 
mal lattice embedded networks, triangles for hollow lattice 
embedded networks with p = 0.5 and circles for p = 1. 
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FIG. 3: (Color online) The average spatial length of connec- 
tions of the normal networks (squares) and hollow networks 
(circles) for scale-free degree distribution (a) and exponential 
degree distribution (b). 
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FIG. 4: (Color online) Random percolation thresholds vs the 
rearrange probability for LESF networks, n = 8. 



FIG. 2: (Color online) Random percolation threshold q c vs 
network side length L for the LESF model (filled symbols) and 
the HLESF model (empty symbols) for different A: Squares 
for A = 2.8 , circles for A = 2.5, and triangles for A = 2.3. 
Each data is the result of averaging 10 4 network realizations. 



increase much compared with the normal networks, for 
small A or large fco- As the average degree (k) becomes 
larger, the difference between (I) decreases. In the limit 
case, the difference goes to 0. While the cost of con- 
structing hollow networks still remains low, they are more 
robust than the normal models. Under the same condi- 
tions, the hollow networks have significant lower perco- 
lation thresholds. 

Based on the above observations, we propose a hol- 
lowing strategy to enhance the robustness of geographi- 
cal networks. For each node in a geographical networks, 
we introduce a probability p to cut down the edges that 
linked to its first n nearest neighbors, then to reconnect 



further nodes in the geographical distance. In this case, 
the degree distribution P(k) deviates a little for small k, 
i.e., around about 10. This only causes a variation in per- 
colation thresholds of a much smaller magnitude. Figure 
^demonstrates the efficiency of the hollowing strategy for 
LESF networks. The percolation threshold drops about 
0.1 with n — 8, namely, it needs 10 percent less nodes 
of the network to maintain a spanning cluster. This can 
have significant effect to prevent the network from break- 
down when the network undergoes a serious crisis and is 
losing its global function. 

For concentrated cases, we study the concentrated 
scale-free networks since the problem mainly concerns 
with immunization strategies [15j |. The long-range links 
(for large A) in the LESF model arc different from those 
in the small- world model [l6j]. In the LESF model, a 
node i first tries to connect to its nearest neighbors; if 
they have already fulfilled their degree quotas, this node 
has to try further nodes until up to the cutoff distance 
A\fki. The long-range links in the LESF model are still 
somewhat local. As far as A is large enough to ensure 
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FIG. 5: (Color online) Properties of concentrated networks, 
(a): Squares are the average degrees, circles are the random 
percolation thresholds, and triangles are the clustering coef- 
ficients, (b): The actual degree distribution for three typical 
A values. The data are log-binned. The assigned degree dis- 
tribution is scale-free, with A = 3.0. 



is indeed the case that, during an epidemic, most people 
are staying home to seclude from the infection. 

To better understand the numerical results, we employ 
the generating function method to determine the per- 
colation threshold for networks with different clustering 
properties. Here the clustering properties can be simply 
depicted by the clustering coefficient, which counts for 
the triangles (3-cycles) in the network, and it can also be 
represented by the number of rectangles (4-cycles) , and 
generally by the number of L-cycles. In the following, a 
general relation of the percolation threshold q c and the 
number of L-cycles for a random network with arbitrary 
degree distribution is determined. As an example, the 
dependence of q c on the clustering coefficient is obtained. 

For uniform occupations (or random failures), the per- 
colation threshold of random tree-like networks is q c — 
(k) I (k(k - 1)) [13, El. It could be obtained by the con- 
dition that the average cluster size diverges, or equiva- 
lently, the average size of clusters that reached by follow- 
ing an edge diverges. A real network is usually clustered 
and contains certain amount of cycles. If the number of 
cycles is small, (e.g. each node belongs to at most one 
cycle) the generating function process can be extended 
to cope with the random percolation problem. 

For a uniform occupation probability q, the generating 
function for the probability of the number of outgoing 
edges of a target node reached by following a randomly 
chosen edge on a clustered network remains the same as 
that of the random tree- like networks [H, IH1 : Fi (%) = 
^■y J2 kP(k)x k ~ 1 . However, if an outgoing edge is not 
independent, i.e., it is terminated to a node that has 
already been visited (having no contribution for reaching 
new nodes) , the node with degree k will only have k — 2 
independent outgoing edges. Thus it is convenient to 
define 
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that almost all the nodes' degree quotas are fulfilled, the 
network properties depend little on A, as shown in Fig. 
[Sfa). However, as A decreases, more and more long-range 
links are prohibited [7]. The network becomes more lo- 
calized, thus the clustering coefficient increases while the 
average degree decreases a little. As A is approximately 
larger than 0.5, most of the nodes' degree quotas are 
still fulfilled except the nodes with large degree quotas. 
Thus the network structure remains mainly unchanged, 
and the percolation thresholds remains almost the same. 
As A decreases further (A < 0.5), a large portion of the 
degree quotas of nodes cannot be fulfilled. Figure [5^b) 
shows a serious degree cutoff (A = 0.3), where the steps 
in the degree distribution come from the symmetry of 
the 2D lattice. Thus the network structure becomes seri- 
ously deteriorated: the average degree decreases sharply, 
the clustering coefficient drops; the percolation threshold 
increases rapidly. These effects indicate that the concen- 
trated networks, are much easier to be immunized. This 



as the generating function of the number of outgoing 
edges reaching new nodes for such a target node. 

Let Hi (x) be the generating function of the size distri- 
bution of the cluster that reached by following an edge, 
by following 2 independent edges, the cluster size distri- 
bution is generated by Hi(x) 2 . But if the 2 edges origi- 
nate from a common node and belong to an L-cycle, the 
generating function should be H[ ll \x)H[ l2 \x), where 
li = \(L — l)/2], I2 = L — 1 — l\. [g] is Gauss' function 

which returns the integer part of g. h[ 1 \x) is the gen- 
erating function for the size distribution of the clusters 
that is reached by an edge and has one edge terminated 
after I steps, and satisfies an iterative relation 

H®(x) = l-F 1 (l)+xF^(H 1 (x))H{ l - 1) (x), 

where the terminal condition H[ (x) is 

h[ 1) (x) = 1-F 1 (1)+xF^ 1 \h 1 (x)). 
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Thus 



H®(x) = (1-^(1)) 



may assume that two triangles could only have at 
most one common node. A node with degree k 
reached by an edge will belong to on average n^{k) = 
(C{k)(k - l)(k - 2) + C{k){k - 1)) /2 = C(k)(k - l) 2 /2 
triangles. Equation Q] will reduce to 



+ (xF[ 1) (H 1 (x))) 1 



In general, we may assume that a node i with degree fcj 
on average belongs to Ui(fc) L-cycles, where nz,(fe) < 1 
or ~ 1. So H%(x) satisfies the self-consistent equation 

H^x) = l-F x (l) + 

^^kP(k)(H[ h \x)H[ h \x)r'<H l (x) k - 1 -^ 

The average cluster size reached by an edge is 
(S) = H[(l) 



(k) 



^2kP(k){n L (k)(Hi h) '(l)+H^'(l)) 



+ (k-l-2n L (k))H[(l)}, 



where 
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fff^l) = <? + «? 
A simple substitution yields 
9 



(S) 



1-9 "W 



1 



^{<fc(* - 1)) - 2 (kn L ) (1 - q<fc g7 a)) )} 



Thus the percolation threshold q c , given by the diver- 
gence of (s), is 



9c = 



(k) 



(k(k - 1)) - 2 (kn L (k)} (l 



2-q l c 1 -q' c 2 q a (k(k-2)) 



2(l-9e) 



<fe> 



(1) 

When 71^ = 0, this result degenerates to the known re- 
sult of the tree-like networks: q c = (k) / (k(k — 1)). In 
general q c < 1, as L — * oo, g' 1 and g' 2 limit to 0, Eq. 
[1] reduces to a second order equation for q c . The root 
with physical meanings is q c — (k) / (k(k — 1)), which is 
just the percolation threshold of the tree-like networks. 
This means that the cycles of infinite length do not affect 
the percolation thresholds. The numerical tests for sev- 
eral typical data values show that, for the same number 
of cycles, q c is an increasing function of g' 1 + q l c 2 . The 
higher the cycle order L is, the less influence it will be. 
Thus the most influential cycles are 3-cycles, which could 
be expressed by the clustering coefficients. 

For L = 3, h = 1% = 1. If the cluster- 
ing coefficient C = (C(k)) is small enough, we 
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{k(k - 1)) - (1 - q c ^^) (C(k)k(k - I) 2 ) 



(2) 



The percolation threshold q c increases monotonically 
with (C(fc)fc(fc— l) 2 ). It is straightforward that when 
C(k) limits to 0, — > 0, q c returns to (k) / (k(k— I)). 
On the other hand, if (C{k)k(k — l) 2 ) diverges, q c maxi- 
mizes to (A;) / (k(k — 2)). Figure [5] shows the dependence 
of the percolation threshold q c on the clustering coeffi- 
cient C. Curves are from Eq. [5] and the symbols are 
from numerical simulation. For simulation, first a ran- 
dom network with prescribed degree distributions is gen- 
erated, then a rewiring process [2l| is applied to achieve 
the required value of clustering coefficient, while keep the 
degree distribution unchanged. Percolation is performed 
using Newman and Ziff's algorithm [131 ]. 
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FIG. 6: (Color online) Percolation thresholds of slightly clus- 
tered networks with truncated degree distribution P(k) ~ 
fc -A e -fc/fc(\ k > m. Squares: A = 3, m = 2, k = 10; Cir- 
cles: A = 6, m = 3, fco — °°- The network size is 10 6 , 
and each data is averaged over 1000 realizations. Solid and 
dashed lines: theoretical result of Eq. ((2)1 for A = 3 and A = 6 
respectively. 

Although Eq. ([2|) holds only for the case of small clus- 
tering coefficient C, the analysis above indicates that for 
large C, or if the network has higher order cycles, the 
fraction of edges that inter-connect existed nodes in a lo- 
cal cluster will further increase. Thus the number of effi- 
cient edges connecting to "new" nodes (in comparison to 
the nodes within the local cluster) may be even smaller, 
which results in a higher percolation threshold. Thus 
when a network is more clustered, or has more cycles — 
not only of order 3 — it will be less robust. 
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FIG. 7: (Color online) Clustering coefficients of the normal 
(squares) and hollow (circles) networks with scale-free (left 
panel) and exponential (right panel) degree distributions. 
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FIG. 8: (Color online) Clustering coefficient spectrum for the 
normal model (squares), hollow model (circles), concentrated 
model with A = 0.5 (up triangles), and concentrated model 
with A = 0.3 (down triangles). The networks are scale-free 
with A = 3.0. The data are log-binned. 



Equation Q] shows that the most influential cycles on 
percolation threshold are triangles; further more, for the 
same number of edges forming cycles, there will be more 
triangles than higher order cycles. Thus the decrease of 
the number of triangles (clustering coefficients) in the 
same kind of networks will be always accompanied with 



the drop of percolation thresholds. Figure [7] displays the 
clustering coefficient C for both normal and hollow net- 
works. It can be seen that the hollow networks have 
much smaller clustering coefficients. Since the percola- 
tion threshold also depend on the spectrum of the clus- 
tering coefficient (Eq. @), we examined the clustering 
coefficient spectrum for the normal, hollow and concen- 
trated scale-free networks, and the results are shown in 
Fig. [8] For normal networks, the spectra for A > 2 are 
almost the same as that for A = 2. The spectra for con- 
centrated networks follow the same scaling law as those 
for normal networks, and the small fluctuations explain 
the behavior of clustering coefficient in Fig. Ufa). For 
hollow networks, the local clustering coefficient is much 
smaller for nodes with small degrees, while it follows the 
same scaling relation as that for normal networks when 
degree is large. Thus the influence of the spectrum of 
clustering coefficients is not crucial. Therefore the de- 
crease in C for hollow networks is consistent with the 
drop of the percolation thresholds (Figs Q] and ^ . 

In short, we have studied how the geographical struc- 
ture affects the percolation behavior of complex net- 
works, and provided analytical understandings by gen- 
erating function formalism on networks with cycles. Our 
study gives a general suggestion on constructing more ro- 
bust real functional networks, such as the Internet, the 
power grid network, etc., that arrange the edges to con- 
nect neighboring nodes as far as possible. Although it 
may cost a little more, it will stand a much reduced risk in 
case of node failures. Also, the hollowing strategy could 
be useful to maintain the global functions of real world 
networks during some emergencies, such as epidemic oc- 
currences, eruptions of electronical virus, or cascade fail- 
ures of power stations, etc. 
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